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Abstract. We examine binary mixtures of superparamagnetic colloidal particles 
confined to a two-dimensional water-air interface both by real-space experiments and 
Monte- Carlo computer simulations at high coupling strength. In the simulations, 
the interaction is modelled as a pairwise dipole-dipole repulsion. While the ratio of 
magnetic dipole moments is fixed, the interaction strength governed by the external 
magnetic field and the relative composition is varied. Excellent agreement between 
simulation and experiment is found for the partial pair distribution functions including 
the fine structure of the neighbour shells at high coupling. Furthermore local crystal 
nuclei in the melt are identified by bond-orientational order parameters and their 
contribution to the pair structure is discussed. 
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1. Introduction 

The mechanisms and principles of heterogeneous crystal nucleation and the subsequent 
microstructure formation are still far from being explored and understood [H El [3]. 
Steering the nucleation behaviour has important implications relevant for protein 
crystallization [4j and the formation of new glasses and metallic alloys 0|6]. Colloidal 
suspensions have been exploited as model systems for crystal nucleation [TllHllQl [TOl [TT] . 
For instance, it is possible to watch crystal nucleation in real space by using confocal 
microscopy [12j. In particular, two-dimensional suspensions of superparamagnetic 
particles confined to the air-water-interface of a pending droplet [13j can easily 
be controlled by an external magnetic field. This allows to tune the interparticle 
interactions and to quench the systems quickly into a supercooled state [Mj. The 
external field induces magnetic dipole moments in the particles p!5] . If the field direction 
is normal to the air- water interface, the resulting dipole moments are parallel and the 
interparticle pair potential is repulsive scaling with the inverse cube of the particle 
separation [16j. Binary mixtures of these superparamagnetic colloids are ideal model 
systems to study crystal nucleation [17j, crystallization [18] and glass formation in 
real-space. 

In this paper, we consider a binary mixture of superparamagnetic particles both by 
real-space microscopy experiment and Monte-Carlo computer simulations of a binary 
dipole-dipole interaction model [20l[l7j. In the simulations, the interaction is modelled 
as a pairwise dipole-dipole repulsion. While the ratio of magnetic dipole moments 
fixed, the interaction strength governed by the external magnetic field and the relative 
composition is varied. We compare the pair correlation functions for strong interactions 
(i.e. for large external magnetic fields) and find good agreement between experiment and 
simulations. Moreover, we discuss the occurrence of peaks in the distance-resolved pair 
correlations in conjunction with local crystallites. The latter are patches of triangular 
and square ordering which are crystalline "seeds" in the amorphous fiuid and building 
blocks from the globally stable crystalline structure [18]. These crystallites actually 
could act as nucleation centers for homogeneous nucleation if the system is quenched 
deeply into the supercooled state [M]. Though heterogeneous nucleation is not tackled 
in the present paper, an understanding of the local crystallites in the bulk is a first 
necessary step in order to access possible pathways of inhomogeneous systems leading 
to heterogeneous nucleation. 

As regards to previous work, we here address strong couplings diflFerent from the 
weakly- coupled case where partial clustering of the small particles was found [20j and 
equilibrated systems. In Ref. [H], a quench was performed and the pair correlation 
functions were found to be widely different from their equilibration counterparts. For 
this genuinely non-equilibrium phenomenon, Brownian dynamics computer simulations 
were employed. 

The paper is organized as follows: in section II briefiy describe the experimental 
set-up and the simulation method. In section III, results of experiment and computer 
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simulations are presented and discussed. Finally we conclude in section IV. 
2. Methods 

2.1. Experimental system and techniques 

The experimental system consists of a suspension of two kinds of spherical and super- 
paramagnetic colloidal particles. The two species are called A and B respectively with 
B referring to the smaller particles. The hard-core diameters of the two species are rf^ = 
4.5 /xm, and rf^ = 2.8 /xm and their magnetic susceptibilities are xa = 6.2 • 10~^^ Am^/T 
and Xb = 6.6 • 10~^^ Am^/T. The details of the experimental set-up are explained 
elsewhere ^71 [I3l |21]. Due to their high mass density, the particles are confined by 
gravity to a fiat water-air interface formed by a pending water drop. The droplet is 
suspended by surface tension in a top sealed cylindrical hole with a diameter of 6 mm 
and a depth of 1 mm of a glass plate. A coil produces a magnetic field H perpendicular 
to the water- air interface which induces a magnetic moment (i.e., = XiH with 
i = A^ B) in each particle. This leads to a repulsive dipole-dipole pair interaction 
[22]. By microscopy, trajectories of all particles in the field of view can be recorded 
over several days providing sufiicient phase space information. The ensemble can be 
considered as ideally two-dimensional as the thermally activated 'out of plane' motion 
of the particles is in the range of a few nanometer. 

While temperature is fixed to room temperature the strength of the interparticle 
interactions is tunable by the external magnetic field strength. A second parameter 
which is varied is the relative composition or the mixing ratio of the particles, 

Nn 

X = — — — (1) 

Na + Nb ^ ^ 

2.2. Monte Carlo simulation technique 

In our Monte Carlo computer simulations, we model the system in two spatial dimensions 
by a pairwise additive potential 

-^^•W = £^ (^,j- = Ai?), (2) 

where r denotes the distance between two particles. For this inverse power potential, at 
fixed composition X, all static quantities depend solely on a dimensionless interaction 
strength - or coupling constant - 

where fc^T is the thermal energy and a = \j s/Pa the average interparticle separation 
between A-particles [23] . Hence effective temperature corresponds to the inverse of the 
coupling r, and the system is completely characterized by three parameters: 1) dipolar 
moment (or susceptibility) ratio 

_ Xb 
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Table 1. Theoretically predicted stable phases [18] for m = 0.1 at T = 0. The same 
notation as in Ref. [18j is used here. The disks (open circles) correspond to A (B) 
particles. 



Phase 



T{A) 



n{A)A3B 



R.e{A)A2B 



K{A)AB 



n{A)A2B2 



S{AB) 



Composition 
(X) 



1/5 



1/4 



1/3 



2/5 



1/2 



Crystalline 
structures 








2) the relative composition X, and 3) the interaction strength F. While we fix the 
former to m = 0.1, the coupling F and the relative composition X are varied. Standard 
Monte-Carlo simulations were performed with A^^ = 400 74-particles and a corresponding 
number Nb of B particles determined by the prescribed relative composition X. The 
particles are in a square box with periodic boundary conditions in both directions. 
Typically 4 x 10^ Monte Carlo steps per particle are used for equilibration and statistics 
is gathered over additional 10^ Monte Carlo steps. 

3. Results 

We present our results for the two considered compositions X = 0.29 and X = 0.44. 
Various coupling strengths F are then investigated by (i) microstructural analysis and 
(ii) partial pair distribution functions. In both cases, real-space experiments and Monte 
Carlo computer simulations have been performed. As a reference, we have gathered 
on table [H the different stable crystalline structures at m = 0.1 from the ground state 
(T = 0) theoretical study [18j. 

3.1. Microstructural analysis 

A visual overview of typical microstructures from the experiments are provided in figure 
[H From the theoretical study [18j (see also table [1]), it is known that the relevant 
stable ground-state crystals consist of pure A-triangular [T(A)] structures {X = 0) and 
intersecting squares of A and B particles [S{AB) phase] at X = 0.5. Local crystallites 
in the fiuid which possess this order are detected by coloring particles which have a pure 
triangular and square order lj| 

:|: In detail, we have used criteria to define A particles which have a pure triangular surrounding of 
other A particles, i.e. which are close to a cut-out of a pure triangular A crystal, and, likewise, we have 
identified A and B particles which form locally an equimolar square lattice S{AB). The corresponding 
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Figure 1. Experimental snapshots for the parameter combination (a) X = 0.29, 
r = 4.9 (b) X = 0.29, r = 38.9 (c) X = 0.29, T = 82.9 (d) X = 0.44, T = 22.6 (e) 
X = 0.44, r = 49.5 (f) X = 0.44, T = 93.9. Big particles are shown in blue if they 
belong to a triangular surrounding and in red if the belong to a square surrounding. 
All other big particles are shown in white color. Few big particles belonging to both 
triangular and square surroundings are shown in pink color. The small particles are 
shown in green if they belong to a square center of big particles, otherwise they appear 
in yellow. 

One can clearly see from figure [1] that the triangular T{A) and square phases S{AB) 
are indeed predominant at strong enough couplings. More precisely, at X = 0.29, see 
figure [H (a-c), there is a strong presence of triangular T(A)-crystallites. This fraction of 
triangular crystallites is growing with increasing F see figure [1] (a-c). 

At nearly equimolarity with X = 0.44, see figure [H (d-f), the situation differs 
qualitatively (compare with figure [H (a-c)), where we have now a strong fraction 

two structure elements are shown in Table [TJ In detail, we associate a triangular surrounding to 
an A particle if the following two criteria are fulfilled simultaneously [17j: 1) The 6-fold bond order 
parameter pq = y^^g^e (where ^6 = | XIatat ^^P (^^^a^a^) with Onn denoting the angles of the six 
nearest neighbour bonds relative to a fixed reference) is larger than 0.94. 2) The relative bond length 
deviation ^6 = | Y1%n where I is the average length of the six bond lengths / nn is smaller than 

0.04. This double condition selects local configurations close to those of a perfect triangular lattice 
where pq is unity and ^6 vanishes. Likewise we define a square surrounding around a B particle by 
the criteria: 1) The 4-fold bond order parameter p^ = ^/^f^ (where ^4 = ^ XIatat^^P (i^^NN) with 
Onn denoting the bond angles of the four nearest neighbour AB bonds) is larger than 0.92. 2) The 
corresponding relative AB bond length deviation 64 is smaller than 0.05. 
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Figure 2. Simulation snapshots for the parameter combination (a) X = 0.29, T = 4.9 
(b) X = 0.29, r = 38.9 (c) X = 0.29, T = 82.9 (d) X = 0.44, T = 22.6 (e) X = 0.44, 
r = 49.5 (f) X = 0.44, r = 93.9. Big particles are shown in blue if they belong to a 
triangular surrounding and in red if the belong to a square surrounding. All other big 
particles are shown in grey color. Few big particles belonging to both triangular and 
square surroundings are shown in pink color. The small particles are shown in green 
if they belong to a square center of big particles, otherwise they appear in grey. 



of squared S(A5)-crystallites. This fraction is increasing with growing F, which is 
consistent with the zero temperature limit predicting the stability of the squared 8(^45)- 
lattice at X = 0.5 (see table [Hand Ref. [18j). 

The simulation snapshots are presented in figure [2] for the same (r,X) parameters 
as in figure [H In a general manner, there is an excellent qualitative agreement between 
the experimental and simulational microstructures, compare figure [1] with figure [2l Very 
interestingly, the theoretically predicted intermediate rectangular phase Ile{A)A2B^ see 
table [H for X = 0.25, is remarkably well present in the snapshot of figure [2](c). This 
feature was not detected in the experiments, compare with figure [I](c), possibly due to 
a slightly imperfect equilibration thereby. 

The emergence of crystalline clusters in the strongly interacting system gives some 
insight into the nucleation behaviour. Presumably the system is in a stable crystalline 
phase at high F but since it is undercooled it does not find the ultimate stable state [24] . 
The intermittent crystal nucleation "self poisons" [25j further crystal growths which may 
be similar to nucleation in liquid crystalline systems [25] , 
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Figure 3. Partial radial pair-distribution functions ^aa(^), gBB{r) and gABi^)- 
Experimental data (EXP) are compared to simulation results (MC) for (a) F = 4.9, 
(b) F = 38.9 and (c) F = 82.9. The composition X = 0.29 is fixed. 



3.2. Pair distribution functions 

We now discuss more quantitatively the structural aspects by inspecting the radially 
averaged partial pair distribution functions, whose corresponding microstructures can 
be found in figured] and figure [2] for the experimental and simulational data, respectively. 

The case X = 0.29 is reported on figure [3l In a general fashion, there is good 
quantitative agreement between experiment and simulation, see figure [3l The only 
situation that slightly deviates from this quality of agreement, concerns the partial pair 
distribution at F = 82.9, see figure [3](c). We explain this with a small drift of 

the colloidal system in the field of view. The whole systems contains more than 100000 
particles and is susceptible to perturbations. The drift induced some shear which may 
constrained here the ordering of the system at F = 82.9 and X = 0.29. The primary 
peaks found at r/a = 1 in S'aaI^) and at r/a = 1/V2 ^ 0.71 in gAuir) is the signature 
of the squared S(A5)-crystallites, see figure [3](b,c). 

The disordered aspect of the material can be best identified by analyzing at the 
partial pair distribution Figured (b,c) shows that the first peak in gBBi"^) is 

located at r/a 0.45, which is much smaller than the (square) unit lattice parameter, 
which is at r/a = 1. This is reminiscent to the similar trend for the small particles to 
form pairs and/or clusters at lower coupling [20j, as can also be easily observed on the 
microstructures of figures [H and [2l 

The case X = 0.44 is reported on figure [H The agreement is now even better 
than at X = 0.29, becoming perfect for all the partial distribution functions (including 
9BB{r)). 
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Figure 4. Partial radial pair-distribution functions ^aa(^), gBB{r) and gABi^)- 
Experimental data (EXP) are compared to simulation results (MC) for (a) F = 22.6, 
(b) F = 49.5 and (c) F = 93.9. The composition X = 0.44 is fixed. 



4. Conclusion 

In conclusion we have put forward the idea that two-dimensional binary mixtures 
are excellent model systems for crystal nucleation as they are easily supercooled by 
increasing an applied external field and crystallize into a variety of crystal structures. 
As revealed by the excellent agreement in the pair correlation functions, the system can 
be modelled by a simple dipole-dipole interaction potential [22j. The strongly interacting 
fiuid bears some crystallites which were identified and could act as possible nucleation 
centers both for homogeneous and heterogeneous nucleation. 

Future investigations should consider the nucleation at fixed imposed nucleation 
seeds which was studied theoretically for one-component two-dimensional systems [26j . 
Steering the nucleation and growth for binary systems is expected to be much richer 
since there are several competing crystalline structures. Exploring more asymmetries 
in the magnetic moments is possible by exploiting nonlinear saturation effects in the 
magnetic susceptibility at high external magnetic fields. This is another parameter 
which is crucially determining the phase behaviour. Finally, it might be interesting 
to use binary charged suspensions confined between charged glass plates [271 ES] as 
two-dimensional model system for crystal nucleation. The interactions are then well- 
approximated effective Yukawa potentials [29l [30j where the screening length is steered 
by the salinity. Again the ground state crystal structures show a wealth of possible 
crystals as recently revealed by lattice-sum calculations [31]. Therefore a rich scenario 
of crystal nucleation and growth phenomena are expected to occur here as well. Another 
realization of dipolar mixtures in 2d are granular systems [32j which show interesting 2d 
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ordering effects . Finally binary colloid mixtures with added nonadsorbing polymers 
will result in effective attractions and possibly liquid-gas phase separation [34[ [35] . 
The interplay of vitrification or gelation and the fluid-fluid phase separation in two- 
dimensions should be an interesting topic for future research [36j . 
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